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1. Introduction 

The decay of J PC = 1~~ charmonium states, 
like J/ip and tp(2S), to K%K% is a CP violating 
process. However, since CP violations in both the 
K°K° and B°B° systems have been well estab- 
lished 1 , it is of interest to search for them in 
other possible processes. 

Furthermore, it has been shown [2] that the 
K g K s system can be used to test the Einstein- 
Podolsky-Rosen (EPR) paradox versus quantum 
theory The space-like, separated coherent quan- 
tum system may yield some J/ip — > K S K S decays 
if EPR's locality is correct while quantum theory 
forbids this process. 

MARK-III searched for decays of J/ip — » 
K S K S with 2.7 million J/ip events [3]. No signal 
was observed, and the upper limit on the decay 
branching ratio was determined to be B(J/ip — > 
K S K S ) < 5.2 x 1CT 6 at the 90% confidence level 
(C. L.). The same search for ip(2S) decays has 
not been performed before. 

In this Letter, we report on a search for J/tp — > 
KgKg using a sample with 20 times more statis- 
tics than before, and on the first search for 
ip(2S) — > KgKg. The data samples used for the 
analyses are taken with the Beijing Spectrome- 
ter (BESII) detector at the Beijing Spectrometer 
(BEPC) storage ring at a center-of-mass energies 
corresponding to Mj/^ and M^sv The data 
samples contain (57.7±2.7) x 10 6 J/ip events and 
(14±0.7) x 10 6 ip(2S) events, as determined from 
inclusive hadronic decays |4I5| . 

BES is a conventional solenoidal magnet detec- 
tor that is described in detail in Ref. [Sj; BESII 
is the upgraded version of the detector 0. A 12- 
layer vertex chamber (VC) surrounding the beam 
pipe provides trigger information. A forty-layer 
main drift chamber (MDC), located radially out- 
side the VC, provides trajectory and energy loss 
(dE/dx) information for charged tracks over 85% 
of the total solid angle. The momentum resolu- 
tion is <7p/p — 0.017\/l +p 2 (p in GeV/c), and 
the dE/dx resolution for hadron tracks is ~ 8%. 
An array of 48 scintillation counters surrounding 
the MDC measures the time-of-flight (TOF) of 
charged tracks with a resolution of ~ 200 ps for 
hadrons. Radially outside the TOF system is a 12 



radiation length, lead-gas barrel shower counter 
(BSC). This measures the energies of electrons 
and photons over ~ 80% of the total solid angle 
with an energy resolution of oe/E = 22%/^/E (E 
in GeV). Outside of the solenoidal coil, which pro- 
vides a 0.4 Tesla magnetic field over the tracking 
volume, is an iron flux return that is instrumented 
with three double layers of counters that identify 
muons of momentum greater than 0.5 GeV /c. 

2. Monte Carlo 

A Monte Carlo (MC) simulation is used for 
the determination of mass resolutions and detec- 
tion efficiencies. For the signal channels, J/ip — > 
K° S K% and ip(2S) -> K°K° S , the angular distri- 
bution of one Kg is generated with the same dis- 
tribution as for the CP allowed channels J/ip — > 
K° S K Q L and ip(2S) -> K° S K° L , namely as sin 2 9, 
where 9 is the polar angle of the Kg in the lab- 
oratory system. Only K s — » ir + ir~ is generated, 
and the K s is allowed to decay in the detector 
according to its lifetime. The simulation of the 
detector response is done using a Geant3 based 
program with detailed consideration of the detec- 
tor performance (such as dead electronic chan- 
nels). Reasonable agreement between data and 
Monte Carlo simulation has been observed in var- 
ious channels tested, including e + e~ —* (j)e + e~ , 
e + e~ — > (7)/!+^, J/ip — > pp and ip(2S) — > 
tt+it-J/iP,J/iP -> £+£- (£ = e,/i). For this 
study, 10,000 events for each signal channel are 
produced. 

3. Event selection 

For the decay channels of interest, the candi- 
date events are required to satisfy the following 
selection criteria: 

1. The number of charged tracks is required 
to be four with net charge zero. Each track 
should satisfy |cos#| < 0.8, where 9 is the 
polar angle in the MDC, and should have 
a good helix fit so that the error matrix of 
the track fitting is available for secondary 
vertex finding. 

2. The tracks are assumed to be either tt + or 
7r~. The higher momentum positive track 
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and the lower momentum negative track are 
assumed to come from one of the Kg de- 
cays, and the remaining two tracks are as- 
sumed to come from the other K s decay. 
The intersection of the track-pairs near the 
interaction point are determined and are re- 
garded as the Kg vertices. 

3. Each K s candidate is required to have 
a decay length in the transverse plane 
(L xy ) greater than 3 mm and an invari- 
ant mass within twice the mass resolu- 
tion as estimated from Monte Carlo sim- 
ulations (the mass resolution in J/ip de- 
cays is 7.1 MeV/c 2 and in ip(2S) decays 
7.9 MeV/c 2 ). 

For ip(2S) decays, backgrounds from i/j(2S) — > 
tt + tt~J/t/j with J/ip decaying into two charged 
tracks are removed by further requiring the 
mass recoiling against the lower momentum pos- 
itive and negative tracks, m rec , to be less than 
3.0 GeV/c 2 . The comparison of the recoil mass 
distributions of signal and background events 
from Monte Carlo simulation is shown in Fig- 
ure^ It can be seen that this requirement is very 
efficient in removing background, and the loss of 
efficiency is small (~ 4%). 

Another possible background is from ip(2S) — > 
7Xc j [J = 0, 2) with X cJ -» K s Kg. Figure H 
shows KgKg missing momentum distributions 
for the signal and background channels, all sim- 
ulated by Monte Carlo. We make an additional 
requirement p m iss < 0.1 GeV/c for tp(2S) selec- 
tion. This removes 94% of \ci background and al- 
most all the Xco background, while the efficiency 
for ip(2S) -> KgKg is about 92% according to 
Monte Carlo simulation. For J/tp decays, since 
r\ c —> KgKg is forbidden by parity conservation, 
no requirement is applied. 

The momenta of the two Kg tracks for events 
passing the above selection are plotted in Fig- 
ures and for J/ip and ip(2S) decays, respec- 
tively, together with the regions predicted for sig- 
nal events by the Monte Carlo simulation. The 
two circles correspond to 

J (Pko - Po) 2 + (p b K a - Po) 2 < na, 
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Figure 1. Recoil mass distributions of Monte 
Carlo simulated ip(2S) -> K S K S (blank his- 
togram) and ip(2S) -> Tr+n-J/ip, J/tp -> £+£- 
(I = e,n) (shaded histogram). The histograms 
are not normalized. The requirement m rec < 
3.0 GeV/c 2 is applied in the selection of ip(2S) 
events. 



where p a K0 and p b K o are the momenta of the two 

Kg candidates in the event; a is the K s momen- 
tum resolution, which is 28 MeV/c in J /if) decays 
and 34 MeV/c in ip(2S) decays, as determined 
from Monte Carlo simulations; and po is the ex- 
pected Kg momentum in J/tp (1.466 GeV) or 
ip(2S) (1.775 GeV) two-body decays. The inner 
circles correspond to n = 1, and the outer circles 
are for n = 2. It can be seen that there are no 
candidate events in either J/tp or ip(2S) decays 
within the one a circle, and there is only one can- 
didate in each case lying between the one and two 
a circles. Using the 2a circle for our event selec- 
tion, one event is obtained each for J/ip — > K S K S 
and for i/j(2S) — > KgKg. The corresponding ef- 
ficiencies, Emc, are 20.74% and 19.18% for J/i/j 
and ip(2S) decays, respectively, as estimated from 
Monte Carlo samples. 

The backgrounds remaining in both J/ijj and 
4>(2S) decays after the K s requirement are stud- 
ied with Monte Carlo simulations. The main 
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Figure 2. Missing momentum distributions of 
Monte Carlo simulated ip{2S) -> K a s K s (er- 
ror bars) and tp(2S) -> -fXcj, XcJ -> K a s K s 
(J = 0,2) (histogram). The histograms are 
not normalized. An additional requirement of 
Pmiss < 0.1 GeV/c is applied to remove the XcJ 
background. 



Figure 4. Momentum of one Kg versus the other 
for the selected ip(2S) decay candidates (open 
squares). The circles correspond to one and two 
a regions of Monte Carlo ip(2S) — > K S K S events. 
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Figure 3. Momentum of one Kg versus the 
other for the selected J/ ip decay candidates (open 
squares). The circles correspond to one and two 
a regions of Monte Carlo J/ip — > ^5^5 events. 



background channel in J/tp decays is from J ftp — > 
K*(892)°K° + c.c. is estimated to be 1.0 ± 0.4 
events, which is obtained using a Monte Carlo 
sample six times as large as the J ftp sample. The 
main backgrounds channels in ip(2S) decays are 
f(2S) -> 1X c2, Xc2 -» K° S K° S and V(25) - 
X*(892)°i ; i: + c.c. A Monte Carlo simulation 
estimates the total background to be 0.67 ± 0.25 
events, with about 80% due to \c2 -> K^K%. The 
errors on the estimated backgrounds are from the 
uncertainties of the branching ratios used and the 
limited statistics of the Monte Carlo samples. 

The two candidate events are investigated fur- 
ther, and it is found that the J/tp — » KgKg can- 
didate is most likely a J/ip -> K~*{892)°K° + c.c, 



with K 



K 7r + misidentified as K9 



because no particle identification is used in the 
event selection, while the ip(2S) — > KgKg can- 
didate looks like a ip(2S) — > 7Xc2, Xc2 — > KgKg 
event with an isolated photon with BSC energy 
of 126 MeV. 
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4. Trigger efficiency 

Due to the long decay length of the high mo- 
mentum K s and the requirement of hits in the 
VC, the trigger efficiency of K%Kg events is lower 
than for normal hadronic events. Since the trig- 
ger system is not included in the Monte Carlo 
simulation, the trigger efficiency is measured with 
data by comparing events within and beyond 
the outer radius of the VC using J/ip(tp(2S)) — > 
KgR^ events, which yields trigger efficiencies of 
(86.7±0.9)% and (81.3 ± 1.9)% for a single K% in 
J/il> and ip{2S) decays, respectively. For KgK s 
events, the trigger efficiency, e tr i g , is calculated to 
be (98.2 ±0.2)% in J/ip decays and (96.5 ±0.7)% 
in ip{2S) decays. 

5. Secondary Vertex Finding 

The secondary vertex finding algorithm is 
checked using J/ip — > K*(892)K + c.c. events. 
It is found that the Monte Carlo simulates data 
fairly well, and extrapolating the difference be- 
tween data and Monte Carlo simulation to the K s 
momentum range under study and correcting by 
the polar angle dependence of the efficiency, cor- 
rection factors of (96.4±3.1)% and (98.1±4.0)% 
are obtained for a single K s in J / tp and ip(2S) de- 
cays, respectively For K S K S decays, the correc- 
tion factors and errors to the Monte Carlo efficien- 
cies, e 2n d, for J/ip -> K S K S and ip(2S) -► K S K% 
are (92.9 ±4.5)% and (96.2 ±5.8)%, respectively. 

6. Systematic error 

The total systematic error on the branching ra- 
tio measurement comes from all sources listed in 
Table, n The simulation of the tracking efficiency 
agrees with data within 1-2% for each charged 
track as measured using channels like J/ip — > A A 
and ip(2S) -> n+ir-J/ip, J/ip -> and 8% 

is taken conservatively as the systematic error for 
the channel of interest. 

The Monte Carlo simulated Kg mass and mo- 
mentum resolutions agree with those of data 
within statistical uncertainties, as has been 
checked with J/ip and ip(2S) -» K S K° L sam- 
ples .8 9 . The requirement that the Kg mass 
and momentum be within two standard devia- 



Table 1 



Summary of the systematic errors. 



Source 


J/1> (%) 


iP{2S){%) 


MC statistics 


2.1 


2.0 


Trigger efficiency 


0.2 


0.7 


2 nd vertex 


4.8 


6.0 


MDC tracking 


8.0 


8.0 


Resolutions 


5.4 


5.6 


Number of events 


4.7 


5.0 


B{K° S -> 7T+7T-) 


0.8 


0.8 


Sum 


12.0 


12.7 



tions introduces systematic uncertainties at the 
5-6% level, dominated by the statistical preci- 
sions of the comparisons between data and Monte 
Carlo simulation. 

The systematic errors on the p m isa and the re- 
coil mass requirements for the ip(2S) sample de- 
pends on the simulation of the momentum of the 
charged tracks and are already included either in 
the tracking or in the systematic error quoted for 
the Kg mass and momentum resolutions. They 
are not further considered here. 

The systematic errors on the total numbers of 
J/ip and ip(2S) events are taken as 4.7% and 
5.0%, respectively, and are measured using inclu- 
sive hadronic events with four charged hadrons 
in the final state 0] for J/ip and using inclusive 
hadrons for ip(2S) The systematic error on 
the branching ratio used, B(K S — ► n + n~) is ob- 
tained from the Particle Data Group ^01 directly. 
Adding all the systematic errors in quadrature, 
the total systematic errors are 12.0% and 12.7% 
for J/ip and ip{2S) decays, respectively. 

7. Results 

Since the observed numbers of events agree 
with the expected background levels for the chan- 
nels under study, upper limits are conservatively 
set by not subtracting background and taking the 
events in MC predicted region as signal. With one 
observed event, the upper limit on the number of 
events is 4.74 at the 95% confidence level. 

The upper limits on the branching ratios of J/ip 
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and ip(2S) 
B(R^K° S K° S )< 



Kg Kg are calculated with 



Ijjl { £ MC • 



N R ■ B(K° S 



where is the upper limit of the number of 
observed events, and Nr is the total number of 
resonance R events. Using the numbers listed 
in Table |3 one obtains the upper limits on the 
branching ratios at the 95% C. L.: 

B(J/ip -» K° S K S ) < 1.0 x 10~ 6 , 

B(i/j(2S) -» K° S K° S ) < 4.6 x 10~ 6 , 

where the systematic errors are included by low- 
ering the efficiencies by one standard deviation. 



Table 2 

Numbers used in the calculations of upper limits 
and the final results. 



R 


J/*l> 


<P(2S) 


n obs 


1 


1 




4.74 


4.74 


SMC (%) 


20.74 ±0.41 


19.18 ±0.39 


etr a (%) 


98.2 ±0.2 


96.5 ±0.7 


£2nd (%) 


92.9 ±4.5 


96.2 ±5.8 




57.7 ±2.7 


14.0 ±0.7 


B(K° S -» 7T+7T-) 


0.6860 ±0.0027 




1.0 x 10" B 


4.6 x 10~ H 



8. Summary 



K Q S K° S and 



The CP violating processes J ftp 
ip(2S) -► K° S K% are searched for using the BESII 
58 million J/tp event and 14 million ip(2S) event 
samples. The upper limits on the branching ra- 
tios are determined to be B(J/ip —> KgKg) < 
1.0xl0- 6 and6(i/'(2S') -> K° S K° S ) < 4.6xl0~ 6 at 
the 95% C. L. The former is much more stringent 
than the previous MARK-III measurement [3], 
and the latter is the first search for this channel 
in ip{2S) decays. Current bounds of the produc- 
tion rates still beyond the sensitivity for testing 
the EPR paradox |2J. 
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